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Abstract 

Compact formulas are obtained for the Casimir energy of a rela¬ 
tivistic perfect fluid confined to a L>-dimensional hypercube with von 
Neumann or Dirichlet boundary conditions. The formulas are con¬ 
veniently expressed as a finite sum of the well-known gamma and 
Riemann zeta functions. Emphasis is placed on the mathematical 
technique used to extract the Casimir energy from a T)-dimensional 
infinite sum regularized with an exponential cut-off. Numerical cal¬ 
culations show that initially the Dirichlet energy decreases rapidly in 
magnitude and oscillates in sign, being positive for even D and nega¬ 
tive for odd D. This oscillating pattern stops abruptly at the critical 
dimension of D = 36 after which the energy remains negative and the 
magnitude increases. We show that numerical calculations performed 
with 16-digit precision are inaccurate at higher values of D. 
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I. Introduction 


If a system has boundary conditions, the inhnite vacuum energy is slightly 
altered compared to the the free continuum case; this leads to a force on 
the boundaries called the Casimir force. In 1948, Casimir calculated the 
attractive force between two conducting plane-parallel plates in vacuum due 
to the zero-point fluctuations of the electromagnetic field. There has been 
an enormous amount of theoretical work on the subject since the pioneering 
work of Casimir (for a general review up to 1997 we refer the reader to 
P, H). The earliest experiment to test Casimir’s calculation was carried out 
by Sparnaay 0] in 1958. The results were inconclusive due to large systematic 
errors and uncontrollable electrostatic forces leading to a 100% uncertainty in 
the results. In 1997, a landmark experiment P, P] using a torsion pendulum 
improved signihcantly on previous results. The most recent experiments 
using atomic force microscopes P] and high precision capacitance bridges 0 
are now in agreement with theoretical calculations to within 1%, eliminating 
any doubt as to the reality of the Casimir force. 

In this work, we calculate the Casimir energy for phonons in a relativis¬ 
tic perfect fluid confined to a H-dimensional hypercube using the cut-off 
method. The Casimir energy of a scalar field in a rectangular cavity with p 
sides of lengths Oi, 02 , ...Op and D — p sides of characteristic length L » at 


was calculated in [10| using the Epstein zeta function regularization scheme. 


It was shown that Neumann and periodic boundary conditions yield a nega¬ 
tive Casimir energy. Determining the sign for Dirichlet boundary conditions 
turned out to be more complicated and was studied in detail in |^, where 
Epstein zeta function regularization was again employed. In |]^ it was shown 
that in a rectangular cavity with p sides of equal length L and D —p sides of 
length » L, the sign of the Dirichlet energy depends on whether p is even 
or odd. For even values of p, the energy is positive when D is less than a 
critical value Dc and negative when D is above Dc- For odd values of p the 
sign is always negative and no critical dimension exists. It was later shown 
r2| that it is possible for the Dirichlet energy to be positive for odd values 
of p if the sides have unequal lengths. 

One alternative to zeta function regularization is the exponential cut-off. 
The cut-off method was employed in to calculate the Casimir energy of 
scalar helds conhned to parallel plates in higher dimensions. In our article. 


2 


















we apply the cut-off method to a perfect fluid conhned to a D-dimensional 
hypercube. We develop a mathematical technique that enables us to ex¬ 
tract the relevant Casimir term from a D-dimensional inhnite sum: one that 
contains the square root of a sum of D squares modihed by an exponen¬ 
tial cut-off term. This technique makes repeated use of the Euler-Maclaurin 
integration formula and a series expansion for the inhnite sum of modihed 
Bessel functions. For both von Neumann and Dirichlet boundary conditions 
we obtain convenient formulas for the Casimir energy as a function of D. The 
formulas are expressed as a single sum of D terms containing the Riemann 
zeta and gamma functions. Numerical calculations show that the Dirichlet 
energy exhibits a clear oscillating pattern up to D = 35: it is positive for 
even D, negative for odd D and its magnitude decreases rapidly. However, 
this oscillating pattern stops abruptly at the critical dimension of D = 36; for 
D > 36, the sign remains negative and the magnitude increases. In contrast 
to the Dirichlet energy, the Neumann energy is negative for all values of D. 

It is instructive to compare the Casimir calculation of a perfect huid to 
that of the open bosonic string. A string embedded in D spatial dimensions 
supports transverse vibrations in D — 1 orthogonal directions. The boundary 
conditions at the two ends of the string, responsible for the Casimir ehect, is 
independent of the dimension D. Therefore, the number of dimensions does 
not complicate the Casimir calculation: the quantity D — 1 contributes only a 
multiplicative factor (see |3 for details). A fluid conhned to a D-dimensional 
hypercube supports longitudinal vibrations in D orthogonal directions. In 
contrast to the string, it has boundary conditions in all D directions leading 
to a Casimir energy with a non-trivial dependence on D. Simply put, for 
the open bosonic string one needs to calculate a single inhnite sum which 
is multiplied D — 1 times whereas for the huid one needs to calculate a 
D-dimensional inhnite sum. This rehects the fact that the perfect huid is 
described by one scalar held which is a function of D-l-1 spacetime dimensions 
whereas the string is described by D -|- 1 scalar helds each a function of two 
spacetime dimensions. 


II. The Acoustic Modes in a Relativistic Perfect Fluid 

A perfect huid is dehned as having at each point a velocity v such that an 
observer moving with this velocity observes the huid as being isotropic. This 
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occurs when the mean free path between collisions is small compared to the 
wavelength. In a frame where the fluid is at rest at some particular position 
and time the energy-momentum tensor has spherical symmetry and is 
given by JT4 

Tij = = 0 ; = p (1) 

where p is dehned as the proper energy density and P the pressure. In the 
same frame the current four-vector is given by 

AT* = 0 ; iV° = n (2) 

where n is dehned as the particle number density. The motion of the huid is 
governed by conservation of energy-momentum and particle number i.e. 

= 0 ; = 0 . (3) 

Small perturbations from equilibrium {p = po, P = Pq and n = no) lead to 
sound waves with the following scalar equation 

d^pi^x' 


dH 


— v^W^p^x) = 0 


(4) 


where p{x) is a scalar held, v is the speed of the sound waves given by 
and X = (x, f). Consider the huid conhned to a D-dimensional 


V = 


- Po 


hypercube with sides of length L. The von Neumann (N) and Dirichlet (Di) 
boundary conditions at a:* = 0 and x* = L are d^p{x) = 0 and p{x) = 0 
respectively (where i=l,2,...,D). The solution to the wave equation (| 
the von Neumann and Dirichlet boundary conditions are respectively, 


for 


{ni}=0 


a 


{ni} 


Awt I ^ —iijjt 

e + a{ni} e 


D 

n 

2 = 1 


,ni7TX\ 

cos(—-—) + at + h 


and 




t 

{rii} 


e + a{ni} e 


{ni}=l 


D 

n 

2=1 


. ,ni7ix\ 
sm (—-—) 


where a; is given by 


(J = ^ {nl + nl-] -h = tt (3 {nl + nl -h . 

Jj 


(5) 


( 6 ) 


(7) 


The parameter (3 = v/L is dependent on the physical and geometrical prop¬ 
erties of the huid: the pressure Pq, the proper density po and the proper 
length L of the sides of the hypercube. 
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III. Quantization and Casimir Energy 

After imposing equal time commutation relations on the scalar field p{x) i.e. 

[p(x, t), p(x', f)] = (x - x') (8) 

one obtains the well known form for the vacuum energy E = (where 

h = 1). For the D-dimensional perfect fluid in consideration, uj is given by 
( 0 ) and the vacuum energy in the Neumann(N) and Dirichlet(Di) cases are 

„ oo 

E-4 Y. (ni + nl + ---+nlf\ (9) 

The multiple sum corresponds to the vacuum energy of the fluid with bound¬ 
ary conditions and is divergent due to the high-frequency modes. The vac¬ 
uum energy with no boundaries i.e. of the continuum, is given by multiple 
integrals and is also divergent. It is the difference between these two energies 
that is of interest and leads to the hnite quantity we call the Casimir Energy 
(the energy needed to set up the boundaries starting from the continuum). 
To extract the relevant constant from the inhnite sum (|^), one regularizes 
the sum to isolate the inhnite contribution of the continuum from the hnite 
contribution stemming from the boundary conditions. There are many ways 
to regularize a sum. In this article we choose an exponential cut-oh term 
g-a(n^+n 2 H hn^} where the parameter a is a positive real number. The 
regularized vacuum energy in the Dirichlet case is then 



In the Neumann case the sums start at Uj = 0 instead of n* = 1. The 
regularized vacuum energy m is hnite and is a function of the parameter 
a. Our goal is to obtain the leading terms in this sum as a —0 and extract 
the constant Casimir term as a function of the dimension D. To accomplish 
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this task we make repeated use of the Euler-Maclaurin integration formula: 
a formula that relates an infinite sum of a function to its integral i.e. 


E/w 


2 = 1 


f{x) dx- - /(O) 




where ^^0) are odd derivatives of / evaluated at zero. There are D 
sums in to evaluate and we apply the Euler-Maclaurin formula to each 
sum except the last one. In Appendix A we show that for the exponential 
function / in (IT), the value of is always zero except for the last 

sum. At the last sum, the value of can diverge and oscillate between 

positive and negative infinity (depending on the value of p) and we therefore 
use a different method of calculation. To summarize, we convert D — 1 sums 


in ([l0|) into multiple integrals by repeated application of the Euler-Maclaurin 
formula and then evaluate separately the last sum. We see from ([TlD that 
each sum (except the last one) gets replaced by an integral of the function 
minus half of the function at zero. This can be expressed by a simple and 
useful prescription 

E-f-l ( 12 ) 


The prescription (IT) can be applied repeatedly to convert multiple sums to 
multiple integrals. The case D = 3 is illustrated below where two of the 
three sums are replaced by (^): 


OO OO OO 

EEE^ 

223 = 1 222 = 1 22i = l 


— a (22^+222+223)^/^ 


E 

713 = 1 


E 

713 = 1 




713 = 1 




^-0.713 


713 = 1 ' 


713 = 1 
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To evaluate (0), we apply D — 1 times the prescription given in (0). This 
yields 


bd. = {i 

nn=l 


nD=- 

p=0 

p=0 


D-1 


= VT/ 


D-l 

p 

D - 1 
P 


= 7Tf3i-l)^2-^da (l-2 f 

nD=l ^ 

C50 pp 

2 )-a. E J 


D-l 


no=l 
-2Y da I{p, a) 


(13) 


where /(p, a) is defined by 
^ rP 

I{p^a) = Y =Zl/ 

n=l n=l 


(14) 


To determine ([T3|) we need to evaluate da I{p, a). The p-dimensional integral 
in I[p, a) can be expressed in terms of the derivative of the Modified Bessel 
function Kp-i (an) B 


e-“("’+-'>+ ...dx^ = -2^ ir'^d^ ( K^{an) (- 


Using the identity 


d 


z dz 


{z-^KYZ)} = {-ir K,+^{Z) 


with z/ = 0, m = Z = an yields 


(-1)^ (] Kn(an) = Kp-i (an) 
' a da J 2 ' 


s 

n\ — 


By substituting (P!7D and ([T3| ) into (EP one obtains 


a,J(p,a)=2^7r"^(-l)^(a„)^ 


d 


a da 


P-1 


s P~1 

n\ — 
a 


Y^oian) 

n=l 


(16) 

(16) 


(17) 


(18) 
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We are interested in obtaining a series expansion of (|^ and isolating the 
relevant constant from the inhnite continuum in the limit as a —0. We 
therefore replace the inhnite sum of Ko{an) by the following series expansion 

Pi 


CXJ ^ 

^Ko{an) = - (C + ln(a/4 7r)) + 

n=l 


2 a 


IT 


E 


I 2 m 71 


(19) 


Consider the terms ln(a/4 7r)/2 and 7r/(2a) in (p^Qf ). They yield terms pro¬ 
portional to 1/aP^^ and respectively in the series expansion of (1^). 

These two terms correspond to the inhnite continuum as a ^ 0. The relevant 
constant related to the Casimir energy stems from the inhnite sum in (|I^ 
i.e. 


lim {da) 

a^O 


d 


a da 


P-1 


E 


71 






TT 


a-0 r(£^)2'^ 


(a^ -|- 7r2)^2 


0{a) 


(-i)^r(^) - 1 


TT 2 2 2 


m=l 

(_i)-^ir(H±2)c(p + 2) 


2p+3 „P+3 
TT 2 2 2 


( 20 ) 


Inserting 


into (|I8D one obtains 
dal(p, a) ^ 


r(£a)C(p + 2) 

2i'+'7r“ 


( 21 ) 


We hnally obtain the Casimir energy for the Dirichlet case by substituting 

m into (0 


D-l 


Em = fd 2(-^-i) (-1)P+^ TT^ r(£^) C(p + 2). 

p=0 


(22) 






















Equation (P^ ) is our final formula for the Casimir energy of a relativistic 
perfect fluid confined to a hypercube with Dirichlet boundary conditions. It 
is conveniently expressed as a hnite sum of D terms involving the gamma and 
Riemann zeta functions; this makes it well-suited for numerical calculations. 
The parameter (3 encompasses the physical and geometrical properties of 
the relativistic perfect fluid: its proper energy density, pressure and length 
L. It plays the same role for the fluid as the string tension does for the 
bosonic string; both (3 and the string tension appear in the Casimir energy 
as dimensionful free parameters. 

Having solved the Dirichlet case it is now relatively straightforward to 
obtain the Neumann case. The regularized vacuum sum in the Neumann 
case, labelled has its sums starting at = 0 instead of Uj = 1 i.e. 

(-a.) E ■ ■ ■ E (23) 

' )^D=0 ni=0 

The above D-dimensional sum can be expressed as a series of fc-dimensional 
sums that start at rij = 1 instead of rij = 0 i.e. fc-dimensional Dirichlet sums. 
The procedure is as follows: we choose k out of the D sums and let these 
k sums start at 1 instead of zero (while the remaining D — k variables are 
not summed and set to zero). One is left with a /c-dimensional Dirichlet sum 
There are (^) ways to choose k among D sums so that the Neumann 
Casimir energy is given by 


D 


(b 


(fc) 

Di 


(24) 


k=l 


where is the fc-dimensional Dirichlet Casimir energy obtained by replac¬ 
ing D by A: in (^^ . Equations (|^) and (|^) are our hnal expressions for the 
Dirichlet and Neumann Casimir energies respectively. 

In table 1 we quote values of the Dirichlet and Neumann Casimir energies 
for D up to 6 calculated using (p2D and (p^). Note that the Neumann Casimir 
energy is negative. In Appendix B, we prove that it is negative for all values of 
D. The Neumann energy, plotted in hgure 1, has a magnitude which increases 
with D. The Dirichlet case is considerably more complicated. Table 1 shows 
that the sign of the Dirichlet energy is negative for odd values of D, positive 
for even values of D and that its magnitude decreases rapidly. These features 
of the Dirichlet energy are valid for low values of D and are plotted in hgure 
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2. The values quoted in table 1 are in agreement with those calculated using 
the Epstein zeta function in (values in [T^ are quoted up to D = 5). It is 
important that the numerical values agree because the Casimir force should 
be independent of the regularization scheme employed. 

Calculations of the Dirichlet energy at higher values of D reveal that 
the oscillation of the sign and the rapid decrease in magnitude stops at the 
critical dimension of D = 36. The Dirichlet energy decreases by twelve orders 
of magnitude from D = 1 to D = 36. To view a plot over such a large span 
requires the energy to be scaled. The magnitude of the Dirichlet energy 
E is less than 1 for the range we consider so that the function |g| io^|e|) 
is well-suited for plotting; it preserves the sign and scales the magnitude 
appropriately. A plot of this function up to D = 110 is shown in hgure 3. 
The distinctive features of hgure 3 are the oscillating pattern which stops 
abruptly at the critical dimension of D = 36 and the ‘plateau’ region which 
emerges immediately afterwards. For D > 36, the energy remains negative 
and the magnitude increases, though slowly in the ‘plateau’ region extending 
to approximately D = 80. To obtain accurate values of the Casimir energy 
at higher values of D, numerical calculations must be performed with greater 
precision than 16-digit precision. We quote in table 2 the Dirichlet energy 
from D = 10 to D = 80 for calculations performed using 16-digit, 24-digit 
and 50-digit precision. With 16-digit precision, numbers begin to show errors 
in the hrst signihcant digit at D = 42 and the sign is wrong for the hrst time 
at D = 49 (yields a positive instead of a negatve sign). Note that in 16- 
digit precision oscillations in the sign resume in the region D > 49. This is 
incorrect; higher-precision calculations show that the sign remains negative 
starting at D = 35. The plot in hgure 3 corresponds to numerical calculations 
done with 50-digit precision; for our plot up to D = 110 this is more than 
enough precision. Note that the 24-digit and 50-digit precision calculations 
yield identical results for values quoted up to D = 80 with four signihcant 
digits. 

We now summarize our results in light of previous work on the Casimir 
energy of scalar helds conhned to rectangular boundaries. One of our goals 
was mathematical: to develop a procedure for calculating (pT]) , a multi¬ 
dimensional inhnite sum regularized with an exponential cut-oh. By repeated 
use of the Euler-Maclaurin formula and a series expansion for the inhnite 
sum of the modihed Bessel function Ko{an), we were able to isolate the 
divergent terms and extract the hnite Casimir energy. In ehect, we reduced 
®) to a hnite sum containing only the gamma and Riemann zeta functions 
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i.e. formula Numerical calculations show that D = 36 is a critical 

dimension, being the hrst even dimension with negative Dirichlet energy. 
Our work focused on the simple geometry of the hypercube whereas previous 
work MM |T^ considered the more general rectangular case and employed 
Epstein zeta function regularization. Results for the rectangular case are 
expressed in terms of asymptotic formulae. In [^, the Dirichlet energy 
for a D-dimensional rectangle with p equal sides is conveniently expressed 
by a single integral with limits running from zero to inhnity and integrand 
containing the elliptic 6 function. In |1l0|, numerical values are quoted 


for low values of D for the hypercube case and they are in agreement with 

for 
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our values. However, the hnite formula (^) was not derived in 
the special case when all sides of the rectangle are equal. 

It is worth noting that the original sum (P) is mathematically a special 
case of a more general class of multiple sums involving arbitrary exponents 

i.e. 


M(s; ai,..., od; ai, • • •, od; c) = ^ {ain°^ + ... + + c) \ (25) 


ni,...,n£)=l 


Using zeta function regularization, Elizalde |Ig obtained explicit formulae 
for (^) expressed as an asymptotic expansion containing the Riemann and 
Hurwitz zeta functions. In contrast to the exponential cut-off method, zeta 
function regularization does not require the introduction of new terms like 
exponentials for convergence; one starts with a convergent sum like (^5|) 
valid for Re s > 0 big enough and then one makes an analytical (usually 
meromorphic) continuation to other values of s. For a detailed mathematical 
treatment of the zeta function regularization theorem and its applications to 
the Casimir energy, the reader is referred to 
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A Appendix 


In applying the Euler-Maclaurin formula ([TI|) to the sums in (p!0D , we show 
that = 0 (except when applied to the last sum in ([lOD ). The most 

general form for the function / is a fc-dimensional integral 




= / Gdx^...dxk (26) 


where G = In f there are k continuous variables 

Xi,... ,Xk which run from zero to inhnity and there are q discrete variables 
rii,... ,nq which run from one to inhnity i.e. / is being summed q times. 
Our goal is to show that the odd derivatives of G with respect to one of the 
discrete variables, say rii, evaluated at ni = 0 is zero i.e. that = 0. 

The hrst derivative of G with respect to ni is G' = —arii G {n\ + ... + + 

xl + ... + xl)~^^‘^ = —ariiGH where H = {nl + .. . + n‘^ + xl + ... +xl)~^^'^. 
The derivative of H with respect to ni is H' = —ni H^. Note that G' and 
H' are expressed in terms of G, H, a and rii. Any subsequent derivatives of 
G' will therefore contain terms of the form 


dn\H^G (27) 

where i, j and I are non-negative integers. Every additional derivative of 
G either increases or decreases j by one. Two consecutive derivatives will 
therefore produce an even change in j. The hrst derivative of G, G', has 
j = 1 so that an additional even number of derivatives applied to G' leads 
to j being odd and positive. Therefore, odd derivatives of G cannot produce 
terms with j = 0. As long as G and do not diverge at rii =0, the terms 
(1^) are zero at rii = 0. Clearly, G does not diverge at rii = 0. does not 
diverge at rii = 0 as long as q > 2 i.e. after ni is set to zero the denominator 
in is never zero if there exists at least one other discrete variable besides 
ni. Therefore, if g > 2, the odd derivatives of G evaluated at ni = 0 are zero 
(and hence the odd derivatives of / evaluated at ui = 0 are zero). 

If Ui is the last discrete variable i.e. g = 1, then if™ evaluated at ni = 0 
diverges at the point where the limits of integration are zero. The last sum 
is therefore calculated using a diherent method. 
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B Appendix 


In this Appendix we show that the Neumann Casimir energy is negative for 
all values of D. The Neumann Casimir energy is given by i.e. 


En — 



E 


(k) 

Di 


(28) 


(k) 

where is the fc-dimensional Dirichlet Casimir energy obtained by replac¬ 
ing D by /c in (^2]) i.e. 




i=l 


k-i\ (-2)*-! r(z)C(* + i) 


i — 1 


22*-i7r^r(f) 


(29) 


Substituting into (psD one obtains 

D k 


En — 


k=l i=l 
D 


D\ fk-1 


k 


i — 1 


- 1 ) 


k+i—l rs-i-k ^ 


vr 2 ((* + 1) 


m 

m 




2=1 


k=i 


-1 


k—i 

30) 


where we used the equality (—= —(—l)^“h Note the change in the 
limits of the double sum i.e. k runs now from i io D and i runs now from 1 
to D. This change does not affect the double sum because one obtains the 
same pairs {k,i). To show that En is negative all we need to show is that 
the sum over k in (|3(]|) is positive. The sum over k is 



where y = x + 1. Clearly derivatives of the polynomial -|- y^~‘^ -l- 
■ ■ ■ -f 1 evaluated at ?/ = 1/2 are positive. We therefore have shown that 
the Neumann Casimir energy is negative for all values of D (assuming the 
dimension D are positive integers). 
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Table 1; Dirichlet and Neumann Casimir energies in units of (3 



D=1 

D=2 

D=3 

D=4 

D=5 

D=6 

Eoi 

-0.131 

0.0415 

-0.0157 

0.00625 

-0.00261 

0.00112 

En 

-0.131 

-0.220 

-0.284 

-0.331 

-0.367 

-0.396 
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Table 2: Dirichlet Energy (in units of (3) for Different Precision Calculations 


D 

16-digit precision 

24-digit precision 

50-digit precision 

10 

4.438e-05 

4.438e-05 

4.438e-05 

11 

-2.035e-05 

-2.035e-05 

-2.035e-05 

12 

9.389e-06 

9.389e-06 

9.389e-06 

13 

-4.360e-06 

-4.360e-06 

-4,360e-06 

14 

2.034e-06 

2.034e-06 

2.034e-06 

15 

-9.538e-07 

-9.538e-07 

-9,538e-07 

16 

4.487e-07 

4.487e-07 

4.487e-07 

17 

-2.120e-07 

-2.120e-07 

-2,120e-07 

18 

1.004e-07 

1.004e-07 

1.004e-07 

19 

-4.770e-08 

-4.770e-08 

-4,770e-08 

20 

2.269e-08 

2,269e-08 

2.269e-08 

21 

-1.085e-08 

-1.085e-08 

-l,085e-08 

22 

5.177e-09 

5.177e-09 

5.177e-09 

23 

-2.488e-09 

-2.488e-09 

-2,488e-09 

24 

1.189e-09 

1.189e-09 

1.189e-09 

25 

-5.754e-10 

-5.754e-10 

-5,754e-10 

26 

2.741e-10 

2,741e-10 

2.741e-10 

27 

-1.345e-10 

-1.345e-10 

-1.345e-10 

28 

6.303e-ll 

6.303e-ll 

6.303e-ll 

29 

-3.205e-ll 
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Figure 1: Neumann Casimir energy as a function of the dimension D 
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Figure 2: Dirichlet Casimir energy at low values of D 
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Figure 3: Scaled Value of the Dirichlet Casimir energy E 
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